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Abstract 

Gluon fusion is the main production mechanism for Higgs particles at the 
LHC. We present the QCD corrections to the fusion cross sections for the Higgs 
boson in the Standard Model, and for the neutral Higgs bosons in the minimal 
supersymmetric extension of the Standard Model. The QCD corrections are 
in general large and they increase the cross sections significantly. In two steps 
preceding the calculation of the production processes, we determine the QCD 
radiative corrections to Higgs decays into two photons and gluons. 



1 Introduction 



The Higgs mechanism is a cornerstone in the electroweak sector of the Standard Model 
[5A^]. The fundamental particles, leptons, quarks and gauge particles, acquire the masses 
through the interaction with a scalar field |jl|]. To accomodate the well-established elec- 
tromagnetic and weak phenomena, this mechanism requires the existence of at least one 
weak-isodoublet scalar field. After absorbing three Goldstone modes to build up the lon- 
gitudinal W^, Zl states, one degree of freedom is left over which corresponds to a scalar 
particle. The properties of the Higgs boson, decay widths and production mechanisms, 
can be predicted if the mass of the particle is fixed . 

Even though the value of the Higgs mass cannot be predicted in the Standard Model, 
constraints can nevertheless be derived from internal consistency conditions • Upper 
bounds on the mass can be set by assuming that the Standard Model can be extended up 
to a scale A before perturbation breaks down and new dynamical phenomena emerge. If 
the Higgs mass is less than 200 GeV, the Standard Model can be extended, with particles 
interacting weakly, up to the GUT scale of order 10^^ GeV, a prerequisite to the renor- 
malization of sin^ from the symmetry value 3/8 down to ~ 0.2 at low energies 
For Higgs masses of more than about 700 GeV, the theory becomes strongly interacting 
already at energy scales in the TeV region . For the large top quark mass found experi- 
mentally p|-p!0|, the requirement of vacuum stability sets a lower limit on the Higgs mass. 
For top masses of 150, 175 and 200 GeV, the lower limits on the Higgs mass are 40, 55 and 
70 GeV, respectively, if the fields of the Standard Model become strongly interacting at 
a scale of about 1 TeV. The lower limits are shifted upwards if the Standard Model with 
weakly interacting fields extends up to energies of the order of the Planck scale. They 
decrease dramatically, however, if the vacuum is only assumed to be metastable 0]. 

The most stringent experimental limit on the Higgs mass in the SAi has been set by 
LEP. A lower bound of 63.9 GeV has been found |Tl| by exploiting the Bjorken process 
Z — > Z*H |12|. The search will be extended to a Higgs mass near 80 to 90 GeV by studying 
the Higgs-strahlung e'^e" — Z* — ZH at LEP2 |T^, T^\. The detailed exploration of the 
Higgs sector in e~^e~ collisions for yet higher masses requires the construction of linear 
colliders 



T5|,|TB|. 



The search for Higgs particles after LEP2 will continue at the pp collider LHC |T7|-p!9|] . 
Several mechanisms contribute to the production of SA4 Higgs bosons in proton collisions 



16]. The dominant mechanism is the gluon fusion process [2C] 



pp gg 



H 



which provides the largest production rate for the entire Higgs mass range of interest. 
For large Higgs masses, the fusion process qq 



2T| becomes compet- 
itive, while for Higgs particles in the intermediate mass range Mz < 



WW, ZZ ^ H ^ 

< 2Mz the 

Higgs-strahlung off top quarks and W, Z gauge bosons |^ are additional important 
production processes. 



Rare decays to two photons will provide the main signature for the search of SAi Higgs 
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particles in the lower part of the intermediate range for masses below about 130 GeV. To 
isolate the narrow 77 signal in the huge 77 continuum background, excellent energy and 
geometric resolution of the 7 detectors is mandatory []18|, [1^. Besides, excellent /i-vertex 



detectors may open the gate to the dominant bb decay mode even though the QCD jet 
background remains very difficult to reject |^5|]. [At the expense of considerably lower rates 
the background rejection can be improved for both reactions by selecting Higgs-strahlung 
events where additional isolated leptons from the associated production of Higgs and top 
or W bosons reduce the QCD background.] Above this mass range, Higgs decays to two 
Z bosons - one Z being virtual in the upper part of the intermediate range - will be used 
to tag the Higgs particle through Z decays into pairs of charged leptons [|18|, |1^. The 
background rejection becomes increasingly simpler when the Higgs mass approaches the 
real-Z decay threshold. At the upper end of the standard Higgs mass range of about 800 
GeV the more frequent decays of the Z bosons into neutrino pairs and jets, as well as 
the WW decays of the Higgs boson, with the W's decaying to leptons and jets, must be 
exploited to compensate for the quickly dropping production cross section. 

Supporting arguments for the supersymmetry extension of the Standard Model are 
rooted in the Higgs sector. Supersymmetric theories provide a natural mechanism for 
retaining light Higgs particles in the background of high GUT energy scales [Q. In the 



minimal supersymmetric extension of the Standard Model [A^iSiSTW] two isodoublet scalar 
fields must be introduced to preserve supersymmetry, leading to two CP-even neutral 
bosons h'^ and a CP-odd neutral boson and a pair of charged Higgs bosons H^. 
The observed value of sin^ 9w has been accurately predicted in this theory [28|, providing 
a strong motivation for detailed studies of this theory ||29|| . 

The mass of the lightest Higgs boson is bounded by the Z mass modulo radiative 



corrections of a few tens of GeV [pO| , [Triviality bounds similar to the SAi Higgs 



sector suggest an upper limit of ~ 150 GeV for supersymmetric theories in general ||32 
The masses of the heavy neutral and charged Higgs particles are expected to be in the 
range between the electroweak symmetry breaking scale and the TeV scale. 

Apart from radiative corrections the structure of the AiSSAi Higgs sector is deter- 
mined by two parameters, one of the Higgs masses, in general m^o, and the angle [3 related 
to the ratio of the vacuum expectation values of the two neutral Higgs fields, tg/5 = V2/V1. 
While the overall strength of the couplings of the Higgs bosons to the SAi particles is 
given by the masses, the mixing angles in the Higgs sector modify the hierarchy of the 
couplings considerably. For example, the coupling of to bottom quarks is strongly 
enhanced for large tg/3 compared with the coupling to the heavier top quarks. Except 
for a small area in the [myio,tg/3] parameter space, Z bosons couple predominantly to 

while the complementary coupling to the heavy Higgs boson is suppressed. The 
pseudoscalar Higgs boson does not couple to the gauge bosons at the Born level. In 
addition, the Higgs particles couple to the SlASy particles, with a strength, however, 
which is essentially set by the gauge couplings. 

The couplings determine the decay modes and therefore the signatures of the Higgs 
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particles. Apart from the small region in the parameter space where the heavy Higgs boson 
decays into a pair of Z bosons, rare 77 and rr decays must be utilized to search for 
the neutral Higgs particles [19[ if b quark decays cannot be separated sufficiently well 



from the QCD background. For large Higgs masses, decays into SUSy particles ^ 
can provide additional experimental opportunities. 

The most important production mechanism for SUSy Higgs particles at hadron col- 
liders is the gluon fusion mechanism, similarly to the SM. Higgs boson production, 

and the Higgs radiation off top and bottom quarks. Higgs radiation off W/Z bosons and 
the WW/ZZ fusion of Higgs bosons play minor roles in the SUSy Higgs sector. 

In the present analysis we have studied in detail the gluon fusion of neutral Higgs 
particles in the Standard Model and its minimal supersymmetric extension. The coupling 
of gluons to Higgs bosons is mediated primarily by heavy top quark loops, and eventually 
bottom quark loops in supersymmetric theories. An extensive literature already exists on 
various aspects of this mechanism. 



The fusion mechanism has been proposed in Ref. for the production of SM. Higgs 



particles at hadron colliders, and has been discussed later in great detail [see Ref. 0, 
p!7| , p!8| , p!9[| for a set of references]. The phenomenological issues for the production of 
Higgs particles in the minimal supersymmetric extension of the Standard Model through 



the gluon fusion mechanism were thoroughly discussed in Refs. All these analyses, 
however, were based on lowest-order calculations. 

Higher-order QCD corrections have first been carried out in Refs. ^ for the 
limit of large loop-quark masses in the Standard Model. Later they were extended to the 
M.SSJ^ Higgs spectrum for this case, however, areas of the parameter space in 



which 6-quark loops are important, are not covered by the approximation. The higher 
order QCD corrections of the fusion cross section for the entire Higgs mass range have 
been given for the Standard Model in Ref. [0] and for its supersymmetric extension in 
Ref. As anticipated, the QCD corrections to the fusion processes are important and 



experimentally significant. Quite generally they are positive and the corresponding K 
factors run up to values of ~ 2. 

Besides the total production cross sections, the QCD corrected transverse momentum 
spectra of the Higgs particles |^ as well as the cross sections for Higgs + jet final states 



, 1^ are of great experimental interest. 



The theoretical analysis of QCD corrections to the gluon fusion of Higgs particles 
involves complicated two-loop calculations; generic Feynman diagrams are depicted in 
Fig. 1. Therefore they have first been performed for the simpler case of Higgs couplings to 
two photons. Fig. 2, for which the virtual QCD corrections are a subset of the corrections 
to the Higgs couplings to gluons. Fig. 3. In the experimentally relevant mass range, the 
QCD corrections to the 77 widths of the SM. and MSSM Higgs bosons are small, of 
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order 55|. In the AiSSAi, special attention must be paid to the kinematical 



range in which the heavy quark-antiquark threshold is nearly mass-degenerate with the 
pseudoscalar state so that non-perturbative resonance effects must be controlled K 



The gluon decay width of the Standard Model Higgs particle has been determined also 
in next-to-leading order; the QCD corrections are positive and numerically important 
36, 47|. The QCD corrections to the rare Higgs boson decay H Z7 [and to the reverse 
process Z H'-f] have been presented in Ref. [HS|; in the mass ranges of experimental 
interest they are tiny, of order a^. [The leading electroweak radiative corrections to the 
Hgg, H'j'y and H'-fZ couplings have been evaluated in the heavy top quark limit to 
O^Gpm'f) [^; they are very small.] 



This paper is divided into two parts. In the first part we will discuss the gluon-gluon 
fusion cross section of the Higgs particle in the SA4 in next-to-leading order QCD. The 
photonic and gluonic partial decay widths of the particle are included in the first part 
of the discussion. The calculations of the production cross section and the decay widths 
have been performed for the entire range of possible Higgs masses. The analytical results 
are summarized in the Appendix in terms of one-dimensional Feynman integrals. In the 
limit where the Higgs mass is either small or large compared to the quark-loop masses, 
the integration can be performed analytically and simple analytical results can be derived 
for the production cross sections and the decay widths. In the second part of the paper 
the analysis will be extended to the CP-even and CP-odd neutral SUSy Higgs bosons. 
To ensure a coherent presentation of the results, some material published earlier in letter 
form will be included in the present comprehensive report. 



2 The Higgs Particle of the Standard Model 
2.1 The Two-Photon Decay Width 

The two-photon decay width of the Higgs boson in the Standard Model, 

is of interest for two reasons. In the lower part of the intermediate mass range of the Higgs 
particle, this rare decay mode provides the signature for the search at hadron colliders 



T9|]. The 77 width determines also the cross section for Higgs production in 77 collisions 



3^. Since the if77 coupling is mediated by triangle loops of all charged particles, the 
precision measurement of the 77 width eventually opens a window to particles with masses 
much heavier than the Higgs mass. If the masses of these particles are generated through 
the Higgs mechanism, the couplings to the Higgs boson grow with the masses, balancing 
the decrease of the triangle amplitude with rising loop mass. As a result, the heavy 
particles do not decouple. However, if the masses of the particles are generated primarily 
by different mechanisms [as in supersymmetric theories, for example], their effect on the 
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i^77 coupling is in general small. 

The decay process ^ 77 proceeds in the Standard Model through W and fermion 
loops, Fig. 2a,b. Denoting the fermionic amplitude hj Af and the W contribution by A^, 
the decay rate is determined by [0, |5l[ 

2 



T{H ^ 77) 



^iV,g5A;(r^)+AH/(rH/) 
/ 



(1) 



where Nc is the color factor, Qf the electric charge of the fermion /. The scaling variables 
are defined by 



m 



H 



and 



m 



Tw 



H 



Aw? 



(2) 



w 



The amplitudes Af and Aw can be expressed as 

Aj{t) = 2[r + (r-l)/(r)]/r2 



-[2r2 + 3r + 3(2r-l)/(r)]/r2 



(3) 



where the function /(r) is given by 

I arcsir^ 

fir 




ITT 



T < 1 
T > 1 



(4) 



If the Higgs mass is smaller than the WW and // pair thresholds, the amplitudes are 
real; above the thresholds they are complex. Fig. 4. Below the thresholds the W amplitude 
is always dominant, falling from (—7) for very light Higgs masses to (-5 - 37rV4) at the 
WW threshold; for large Higgs masses the W amplitude approaches Aw (~2). Quark 
contributions increase from 4/3 for light Higgs masses (compared with the quark mass) to 2 
at the quark-antiquark threshold; far above the fermion threshold, the amplitude vanishes 
linearly in r mod. logarithmic coefficients, Af — [log(4r) — 'i7r]^/2r, i.e. proportional 
to rnj/m\. The contribution of the W loop interferes destructively with the quark loop. 
For Higgs masses of about 600 GeV, the two contributions nearly cancel each other |^ . 

Since the Hff coupling is proportional to the fermion mass, the contribution of light 
fermions is negligible so that in the Standard Model with three families, only the top 
quark and the W gauge boson effectively contribute to the 77 width. Since the W and 
fermion loops interfere destructively, a fourth generation of heavy fermions would reduce 
the size of the if77 coupling. For small Higgs masses the additional contributions of the 
heavy quarks and the charged lepton would suppress the decay width by about one order 
of magnitude. 

To fully exploit the potential of the 77 decay mode of the Higgs particle and the 
production in 77 collisions, the QCD corrections must be shown to be under proper 
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control. To include the gluonic QCD corrections, twelve two-loop diagrams plus the 
associated counter terms must be taken into account. Generic examples are depicted in 
Fig. 2c. 

Throughout this analysis we have adopted the on-shell renormalization scheme which 
is convenient for heavy quarks. In this scheme the quark mass mg is defined as the pole 
of the propagatorQ, related in the following way to the running mass 

-| 12/{33-2Np) 

{l + Oial)} (5) 



a 



aJm% 



at the mass renormalization point /ig. It should be noted that this definition of the running 
mass does not coincide with the running MS mass. The wave function is renormalized 
[^2^^] such that the residue at the pole is equal to unity. The photon-quark vertex 
is renormalized at zero-momentum transfer; the standard QED Ward identity renders 
the corresponding renormalization factor equal to the renormalization factor of the wave 
function. Since the fermion masses are generated in the Standard Model by the interaction 
with the Higgs field, the renormalization factor associated with the Higgs-quark vertex 
[Zhqq] is fixed unambiguously by the renormalization factors for the mass and Z2 for 
the wave function. From the Lagrangian 

jj 

Co = -moQoQo — 

V 

H H 

= -niQQQ— + ZHQQmQQQ— (6) 

we find m 

Zhqq = 1 — ^2^m (7) 

In contrast to the renormalized photon-fermion vertex, the scalar HQQ vertex r{p',p) is 
renormalized at zero momentum transfer by a finite amount 7^ of order after subtract- 
ing Zhqq due to the lack of a corresponding Ward identity. The finite renormalization 
7m, corresponds to the anomalous mass dimension discussed later. 

We have calculated the two-loop amplitudes using dimensional regularization. The 
five-dimensional Feynman parameter integrals of the amplitudes have been reduced an- 
alytically down to one-dimensional integrals over polylogarithms ||5^ which have been 
evaluated numerically^ [see Appendix A]. In the two limits where m'jj/AmQ is either very 
small or very large, the amplitudes could be calculated analytically. 

The QCD corrections of the quark contribution to the two-photon Higgs decay am- 
plitude can be parameterized as 



An = A, 



LO 



IT 



(8) 



"'^We have chosen mt — 174 GeV for the t pole mass ^ and rub — 5 GeV for the b pole mass. 

^The scalar integral associated with the gluon correction to the HQQ vertex has also been analyzed 
by means of analytical | |55[ | and novel approximation methods [ p6[ . The results are in agreement within 
an accuracy of 10^^. 
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The coefficient Ch splits into two parts, 

C^ = Ci + C2log^ (9) 

where the functions q depend only on the scaling variable r = m|^/4mQ(/XQ). The same 
running quark mass mQ{^Q), evaluated at the renormalization scale /zg, enters in the 
lowest -order triangle amplitude Aq^. The scale in is arbitrary to this order; however, 
in practice it should be defined of order rriH- As a typical renormalization scale we have 
chosen /xq = m///2. This choice suggests itself for two reasons. The QQ decay threshold 
is [perturbatively] defined at the correct position 2mQ{mQ) = 2mQ. In addition, it turns 
out a posteriori that all relevant large logarithms are effectively absorbed into the running 
mass for the entire physically interesting range of the scaling variable r. 

The correction factor Ch is displayed in Fig. 5, illustrating the preferred choice /ig = 
for the renormalization scale. The coefficient is real below the quark threshold 
and complex above. Near the threshold, within a margin of a few GeV, the present 
perturbative analysis is not valid. The formation of a P-wave O^"*" resonance, interrupted 



however by the rapid quark decay |53], modifies the amplitude in this range [^]. The 
perturbative analysis may nevertheless account for the resonance effects in a dual way. 
Since QQ pairs cannot form 0"*""^ states at the threshold, Ch vanishes there. ^cCh 
develops a maximum very close to the threshold. 

The QCD-corrected 77 decay width of the Higgs boson is shown in Fig. 6a. The 
correction relative to the lowest order is small in general. Fig. 6b. The corrections are 
seemingly large only in the area where the destructive W- and Q-loop interference makes 
the decay amplitude nearly vanish. 

The Limit of Large Quark-Loop Mass 

In the limit m'jj/AmQ 0, the five-dimensional Feynman parameter integrals can be 
evaluated analytically. The correction factor for the H'j'j coupling 

ml/4ml^0: I + Ch—^1-— (10) 

TT TT 

agrees with the result of the numerical integration in this limit. 

The H'j'y coupling can also be derived by means of a general low-energy theorem for 
amplitudes involving soft Higgs particles |T^, limp^^o A{X H) = {mo / v)dA{X) / drriQ. 
The theorem is easy to prove. For zero 4-momentum the kinetic derivative term in 
the Lagrangian can be neglected and the [space-time independent] Higgs field can be 
incorporated by adding the potential energy to the bare mass term, mg ^«o(l + H/v), 
in the Lagrangian. The expansion of the bare propagators for small values of H/v is then 
equivalent to inserting a zero-momentum Higgs field in an arbitrary amplitude A{X), 

^__^_mo^_ ^^^^ 
^ — niQ ^ — mo V ^ — mo 
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and generating this way the amphtude A{XH). Since the bare mass mg and the renormal- 
ized mass mg are related by the anomalous mass dimension, c? log mo = (1 + 7^)*^ logmQ, 
we find for the final form of the theorem 



lim A(XH) 



rriQ d 



-AiX) 



(12) 



It is well-known that the theorem can be exploited to derive the H'j'j coupling in 



lowest order [|T^, However, the theorem is also valid if radiative QCD corrections 
are taken into account. For large fermion masses, the vacuum polarization of the photon 
propagator at zero momentum is given by 



n 



TTlr 



a.. 



l + -^r(l + e) 



'47r/i^ 



1717 



0{e) 



so that 



dli 

dmo 



2^ 

7r 



1 + — 

TT 



From the anomalous mass dimension to lowest order, 
one readily derives the correction Ch of the if 77 coupling 

a., 1 + Cts/TT 



m^/4m. 







1 + C 



H- 



TT 



TT 



(13) 



(14) 



(15) 



(16) 



1 + 2a,/7r 

Compared with the radiative QCD correction to the photon propagator, (1 + a^/Tr), just 



the sign of the correction is reversed, (1 — a^/Tr), for the H'yy coupling In the 



notation of eq.(^ the correction is attributed to ci while C2 ~ l/'^g vanishes for large 
quark masses. 

The same result can be derived by exploiting well-known results on the anomaly in 
the trace of the energy-momentum tensor [^, 



e 



MM 



;i + 7m)"^oQoQo + ^—Ff^uF^u 
A a 



(17) 



Pa denotes the mixed QED/QCD /? function defined by da{fj,^)/d log fi = Pa- Since 
the matrix element (77 1 6^/^10) vanishes for infrared photons, the coupling of the two- 
photon state to the Higgs source (?tio/v)QoQo is given by /3^/[4a(l + 7^)], with = 
2eQa^/7r(l + as/vr) including only the heavy quark contribution to the QED/QCD /? 
function. Thus the H'j'y coupling is described by the effective Lagrangian 



£(i777) 



end 



271 



1/2 



1 _ ^ 

TT 



F F H 
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which is apparently equivalent to the previous derivation of the ^^77 coupling in the limit 

The Limit of Small Quark-Loop Masses 

In the limit mq^iiQ) the leading and subleading logarithms of the QCD correction 
Ch can be evaluated analytically: 



: 



Ch -T^log'f-4r-2e) 



^log(-4r-«e)+21og^ (19) 
6 rrin 



and, split into real and imaginary parts, 



QuiCh 
^cCh 



n 
3 



-log(4r) + 2 



1^ 



log'(4r) 



7r 



2 

^log(4r) + 21og^ 
3 rrin 



(20) 



The choice of the renormalization scale /ig is crucial for the size of Ch- Choosing the 
on-shell definition fiQ = rriQ leads to very large corrections in the imaginary as well as the 
real part, as demonstrated in Fig. 5. By contrast, for /ig = ^mn, the corrections in the 
real and imaginary part remain small in the entire r range of interest, r < a few x 10^ 
for rrih ~ 3 GeV and m^/ < 1 TeV. [This coincides with the corresponding observation 
for the decay H —y bb where the running of the 6-mass up to the scale ^mn absorbs the 
leading logarithmic coefficients 



Only for log r values above the physical range must 
the leading logarithmic corrections be summed up; such an analysis is beyond the scope 
of the present investigation. 



2.2 The Gluonic Decay Width 

Gluonic decays of the Higgs boson 

H^gg 

are of physical interest for arguments similar to the preceding section. However, there 
are some qualitative differences. Since the particle loops mediating the Hgg coupling 
carry color charges, the color-neutral W, Z gauge bosons do not contribute. The gluonic 
branching ratio can only be measured directly at e^e~ colliders and for Higgs masses 
presumably less than about 140 GeV |l^ since it drops quickly to a level below 10~^ 
for increasing masses. In this range, a fourth generation of fermions would enhance the 
branching ratio to a level where it becomes competitive with the dominant bb decay mode. 

The gluonic width determines the production cross section of Higgs bosons in gluon- 
gluon fusion to leading order at hadron colliders. The cross section, however, is strongly 
affected by QCD radiative corrections so that the width can be measured in this indirect 
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way only within about 20%. 

At the Born level the contribution of heavy quarks to the gluonic width in the Standard 
Model is given by 



m ^ gg) 



GfolI 
367277' 



-m 



H 



(21) 



where Aq denotes the quark amplitude, already discussed in eq.(^, without the color 
factor. The top quark contribution is by far dominant in the SM.. Any additional heavy 
quark from a fourth family etc. increases the decay amplitude by a factor 2 in the limit 
where the Higgs mass is small compared with the QQ threshold energy. 

The QCD corrections to the gluonic decay width ^ are large. Several classes 
of diagrams must be calculated in addition to those familiar from the two-photon decay 
amplitude. Generic examples are shown in Fig. 3. The virtual corrections involve the non- 
abelian three-gluon and four-gluon couplings, and the counter terms associated with the 
renormalization Zg — 1 = {Zi — 1) — |(^3 — 1) of the QCD coupling. We have defined as 
in the MS scheme with five active quark flavors and the heavy top quark decoupled . 
Besides the virtual corrections, three-gluon and gluon plus quark-antiquark final states 
must be taken into account, 

H ggg and gqq 

In the quark channel we will restrict ourselves to the light quark species which we will 
treat as massless particles^]. As a consequence of chirality conservation, the gluon decay 
amplitude does not interfere in this limit with the amplitude in which the qq pair is 
coupled directly to the Higgs boson [g{Hqq) of order rriq, but kept non-zero]. This would 
be different for top quark decays H ttg where the decay mechanism of Fig. 3, however, 
is a higher-order effect, suppressed to 0{g1) already at the amplitude level with respect 
to the gluon bremsstrahlung correction of the basic ti decay amplitude. The light-quark 
final states in the QCD corrections to the gluonic decays, on the other hand, must be 
taken into account since they are energy-degenerate with the gluon final states. 

The result can be written in the form 



^{H gg{g), gqq) = Tlo{H 



a. 



TT 



(22) 



with 



T 



95 7,, 33 -2A^^ , 
— - -Np + — - log Ar + AE 



4 6 



6 



m 



(23) 



H 



■^We include c,b quarks among the light quarks so that all large logarithms log mf^/m'^ ij, associated 
with final state particle splitting, are removed by virtue of the Kinoshita-Lee-Nauenberg theorem. This 
assumes that when the theoretical prediction will be compared with data, c and b quark final states in 
coUinear configurations are not subtracted. Note that in Higgs decays to c, 6 quark pairs plus an additional 
gluon jet, the heavy quarks are emitted preferentially back-to-back and not in coUinear configurations. 
[A more detailed phenomenological analysis of these final states is in progress.] 
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The first three terms survive in the limit of large loop masses while AE vanishes in this 
limit, /i is the renormalization point and defines the scale parameter of ag- It turns out a 
posteriori that the higher order corrections are minimized by choosing the pole mass mg 
for the renormalized quark mass; this is evident from Fig. 7a. The correction AE, given 
explicitely in the Appendix, is displayed in Fig. 7b for the physically relevant mass range. 
In Fig. 8 we present the gluonic width of the Higgs boson including the QCD radiative 
corrections^. 

The total decay width and the branching ratios of all decay processes in the Standard 
Model are shown in Fig. 9 for Higgs boson masses up to 1 TeV. All known QCD and 
leading electroweak radiative corrections are included. 

The size of the QCD radiative corrections depends on the choice of the renormalization 
scale /i for any fixed order of the perturbative expansion. A transparent prescription is 
provided by the BLM scheme in which the Np dependent coefficient of the correction is 
mapped into the coupling as, summing up quark and gluon loops in the gluon propagators. 
We shall apply this prescription in the large loop-mass limit where the amplitude can be 
calculated analytically: 



ml /Ami 



: 



E( 



95 

T 



7 

6^ 



33-2Nf , 

log 



6 



m 



H 



Choosing 

fJ'BLM = e'iniH ~ O.nrriH 
the Np dependent part drops out of E{t) and we are left with 



ggig) + am) = ^eWs^i^'Blm)] 



1 + 



9 

2 TT 



(24) 
(25) 

(26) 



A large fraction of the total QCD correction is thus to be attributed to the renormalization 
of the coupling. 

We shall conclude this subsection with a few comments on the effective Hgg La- 
grangian [^. In the same way as for -^77, we can derive the effective gluon Lagrangian 
for quark-loop masses large compared to the Higgs mass by taking the derivative of the 
gluon propagator with respect to the bare quark mass for = 0. Introducing again the 
anomalous mass dimension 7^, one finds for the gauge-invariant Lagrangian, 



1 



Hgg 



4 l + 7„(a. 



TT 



(27) 



where 



3^ 



1 + 



19 a, 



TT 



and 



7n 



.OLs 



TX 



f28) 



^In all numerical analyses and figures, the contributions of the b quark loops have been included. Even 
for small Higgs masses these effects remain less than about 10% of the leading t quark contributions. 
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so that to second order 



4 TT . 



Gl^G^, H (29) 



As a consequence of the non-abehan gauge invariance the Lagrangian describes, besides 
the Hgg couphng, also the Hggg and Hgggg couphngs, Fig. 10a. 

In contrast to the effective H'yy Lagrangian, Cngg does not describe the Hgg inter- 
action to second order in in totaL This Lagrangian accounts only for the interactions 
mediated by the heavy quarks directly, but it does not include the quantum effects of the 
light fieldsQ: Cngg must be added to the light-quark and gluon part of the basic QCD 
Lagrangian, and this sum then serves as a new effective Lagrangian for Higgs-gluon-light 
quark interactions. Physical observables associated with the low-energy Higgs particle 
are calculated by means of this effective Lagrangian in the standard way, generating gluon 
self-energies, vertex corrections, gluon-by-gluon scattering, gluon splitting to gluon and 
light quark pairs, etc. In summa, the diagrams displayed in Fig. 10b must be evaluated, 
taking into account also the corresponding counter terms that renormalize the coupling 
as and the gluon wave function. 

The fixed-order program discussed so far can be applied to the mass region where 



mH/2mQ is small [in essence < 1] |3^, ^ but logmq/mH still moderate so that 



logarithmic terms need not be summed up. This is the kinematical region of physical 
interest. Based on a careful RG analysis [^, the logarithmic terms have been summed 
up in the limit where also log mg/m/f is large. This leads to the plausible result that 
the energy scale in the effective Hgg Lagrangian is set by the heavy-quark mass while 
the Higgs mass is the scale relevant for the additional light-quantum fluctuations. This 
can be incorporated by substituting asEij) [ll/2]as{mQ) + [73/4 — 7/6NF]as{mH) 
in eg . (plD , leaving us with the light-quantum fluctuations as the main component of 
the QCD corrections in this mathematical limit. For moderate values of logmg /tuh the 
splitting is of higher order in the QCD coupling and can be neglected. 



2.3 Higgs Boson Production in pp Collisions 



Gluon fusion |pO | 



pp-^ gg ^ H 

is the main production mechanism of Higgs bosons in high-energy pp collisions throughout 
the entire Higgs mass range. As discussed before, the gluon coupling to the Higgs boson 
in the Standard Model is mediated by triangular loops of top quarks. The decreasing 
form factor with rising loop mass is counterbalanced by the linear growth of the Higgs 



^Technically, the additional contributions are proportional to a common factor [fi^ /m^^Y which van- 
ishes if the Higgs mass is set to zero before e is driven to (—0). [Note that these mass singularities are 
regularized formally for e < 0.] However, keeping the Higgs mass non-zero but small, the expansion in e 
gives rise to log /i/m/f terms which fix the size of the renormalization scale of the physical process. 
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coupling with the quark mass. [Heavier quarks still, in a fourth family for instance, would 
add the same contribution to the production amplitude if their masses were generated 
through the standard Higgs mechanism.] 

To lowest order the parton cross section. Fig. la, can be expressed by the gluonic 
width of the Higgs boson. 



87r2 



m 



Tlo{H gg) 



(30) 



H 



where s is the gg invariant energy squared. Recalling the lowest-order two-gluon decay 
width of the Higgs boson, we find 



288^271 



(31) 



The tq dependence of the form factor has been given in eq.(|^). With rising mass, the 
width of the SM. Higgs boson quickly becomes broader. This effect can be incorporated in 
the lowest-order approximation by substituting the Breit-Wigner form for the zero-width 
^-distribution 



1 



sTn/mH 



and changing kinematical factors — > s appropriately. 
Denoting the gluon luminosity as 



(32) 



d£9a 
dr 



1 dx 



T X 



g{x,M^)g{T/x,M' 



(33) 



the lowest-order proton-proton cross section is found in the narrow-width approximation 
to be 

o-Loipp ^ H) = cfoTh-^ (34) 



dr. 



H 



where the Drell-Yan variable is defined, as usual, by 



m 



th 



H 



(35) 



with s being the invariant pp collider energy squared. The expression r^dC^^ /drn is only 
mildly divergent for Th ^ 0- 

The QCD corrections to the fusion process gg ^ H |3^, Fig. lb, 
gg H{g) and gq Hq, qq ^ Hg 
14 



involve the virtual corrections for the gg H subprocess and the radiation of gluons in 
the final state; in addition, Higgs bosons can be produced in gluon-quark collisions and 
quark-antiquark annihilation. These subprocesses contribute to the Higgs production at 
the same order of a.^- The virtual corrections modify the lowest-order fusion cross section 
by a coefficient linear in a.^- Gluon radiation leads to 2-parton final states with invariant 
energy s > mf^ in the gg, gq and qq channels. The parton cross sections for the subprocess 
i + j —>■ H + X may thus be written 



a. 



TT 



a. 



(5(l-f) + A,(r,rQ)-0(l 

TT 



for i,j = g, q, q. The new scaling variable f, supplementing the variables th 
tq = m^j/Amq introduced earlier, is defined at the parton level. 



(36) 
m\j/s and 



H 



(37) 



The quark-loop mass is defined as the pole mass in the scaling variable tq. The coef- 
ficients C{tq) and Dij{f,TQ) have been determined by means of the same techniques as 
described for the H'j'j and Hgg couphngs at great detail. The lengthy analytic expres- 
sions for arbitrary Higgs boson and quark masses are given in the Appendix in the form 
of one-dimensional Feynman integrals. The quark-loop mass has been defined in the on- 
shell renormalization scheme, while the QCD coupling is taken in the MS scheme. If all 
the corrections ( p6| ) are added up, ultraviolet and infrared divergences cancel. However 
collinear singularities are left over. These singularities are absorbed into the renormaliza- 



tion of the parton densities We have adopted the MS factorization scheme for the 
renormalization of the parton densities. The final result for the pp cross section can be 
cast into the form 



a{pp ^H + X) = oo 



l + C 



a. 



TT 



Th-: h Aagg + Aagg + A(Tgg 



(38) 



with the renormalization scale in and the factorization scale of the parton densities to 
be fixed properly. 

The coefficient C{tq) denotes the contributions from the virtual two-loop corrections 
regularized by the infrared singular part of the cross section for real gluon emission. 
This coefficient splits into the infrared part vr^, a logarithmic term depending on the 
renormalization scale and a finite rg-dependent piece c{tq), 



C{rQ) = TT^ + c(rQ) 



33 



6 



2Nf /i^ 
-log- 



m 



(39) 



H 



The term c{tq) can be reduced analytically to a one-dimensional Feynman-parameter 
integral [see Appendix B] which has been evaluated numerically . In the heavy-quark 
limit Tq = m'jj/4mQ -C 1 and in the light-quark limit tq ^ 1, the integrals could be 
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solved analytically. 

The (non-singular) hard contributions from gluon radiation in gg scattering, gq scat- 
tering and qq annihilation depend on the renormalization scale /i and the factorization 
scale M of the parton densities [Fig. lb], 



99 



1 , daa a 
dr 



X — i -rPgg{f) log — + dgg{f, Tq) 



TT 



+ 12 



logfl-f) 



1 - f 



91 



Aa. 



99 



•^^H q,q 



dai a, r.p ... 

da^ as , , 

-j^ X —do dgg{T, Tq) 



f[2-f(l -f)]log(l-f 

- dgq{f,TQ 



1 M2 , ^ 

--log — + log(l - r 

2 s 



(40) 



with f = th/t. The renormalization scale enters through the QCD coupling as{fi'^) in 
the radiative corrections and the lowest-order parton cross section cro[as(Ai^)]- Pgg and 
Pgg are the standard Altarelli-Parisi splitting functions pB |, 



P99ir) 



6 



1 



41 + (l-f)^ 
3 f 



;r - 2 + f(l - f) ^ + 



33 - 2Ar^ 
6 



5(1 -f) 



(41) 



F+ denotes the usual + distribution such that F{t)^ = F{f) — 5{1 — f) df'F{f'). The 
coefficients dgg, dgg and dqq can be reduced to one-dimensional integrals [Appendix C] 
which have been evaluated numerically [^ for arbitrary quark masses. They can be 
solved analytically in the heavy and light-quark limits. 

In the heavy-quark limit the coefficients c{tq) and dij{f,TQ) reduce to very simple 
expressions [^, |3^ 



= ^H/^^rnl < 1 : 



c(t-q) 
dgqi'r,TQ) 

dqq{T,TQ) 



11 
Y 



2 ^ ' 



32 
27' 



2r-^- 
3 



(42) 



The corrections of 0{tq) in a systematic Taylor expansion have been shown to be very 
small [^. In fact, the leading term provides an excellent approximation up to the quark 
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threshold ttih ^ '^^mQ. 

For the sake of completeness we quote the differential parton cross sections for hard- 



gluon radiation and quark scattering in the heavy quark-loop limit ^ 



with 



H{gg 
H{gq 



da 
di 



Hg) 
Hq) 



288V27r2 



H{s,t) 



(43) 



3 s^ + i^ + u^ + m^ 



H 



H{qq^Hg) 



2 sHu 
_2 s^ + u" 

3 si 

16 i^ + u^ 
9 s2 



The Mandelstam variables t, u are the momentum transfer squared from the initial partons 
gg, gq, qq, respectively, to the Higgs boson in the final state. [The singularities for £, -u — > 
can be regularized in n dimensions.] 

In the opposite limit where the Higgs mass is very large compared with the top mass, 
a compact analytic result can be derived, too: 



c(i-q) 



5 4 

— log^(-4rQ - ^e) - - log(-4rQ - ie) 



-- log(4TQ) {7 - 7f + 5f2} - 6 log(l - f) {l - f + f2} 

+2i^|3-6f-2f^ + 5f3-6fn 
1 — r ^ J 



d 



qq 







i + (i-f) 



^log(4rQ)+logf^-^ 



(44) 



These approximate expressions are valid to leading and subleading logarithmic accuracy. 

The final results of our analysis are presented in Fig. 11 and the subsequent figures 
for the LHC energy ^/s = 14 TeV. [A brief summary is also given for 10 TeV.] They 
are based on a top-quark mass of 174 GeV [B-|10[- If not stated otherwise, we have 



adopted the GRV parameterizations ||6^ of the parton densities. These are defined in the 
MS scheme^. We have chosen af\mz) = 0.117 of the MS scheme in next-to-leading 

^We may switch to different schemes by adding the appropriate finite shift functions |^ to the 
integrals Auy , i.e. substituting Pijir) logM^/s Pij{T) \ogM^ /s + /y(f) in eqs.(pO[). 
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order. This corresponds to the average measured QCD couphng for five quark degrees of 



freedom with = 214 MeV; the standard matching conditions |^ are adopted 



at /i = : af\fi = mt\A^) = af\^i = mt\k%) with a|L = 0.413 A^. The GRV 
fits are based on a somewhat smaller value of as- This introduces a slight inconsistency 
into the numerical evaluation of the cross section which we allow for since, on the other 
hand, the basic parton cross section is quadratic in as and thus depends strongly on the 
choice of the QCD coupling. In order to correct the difference in the Ajg^ values, the 
factorization scale M at which the parton densities are evaluated, has been changed to 
adjust appropriately the ratio M^/A|j^ which enters in the structure functions^. The 
cross section is sensitive to gluon and quark densities down to x values of order 10~^ to 
10~^, so that subtle non-linear effects in the evolution at small x need not be taken into 
account yet. 

We introduce K factors in the standard way, 

K,.. = ^ (45) 

The cross sections ano in next-to-leading order are normalized to the cross sections aio^ 
evaluated consistently for parton densities and in leading order; the QCD NLO and 
LO couplings are taken from the GRV parameterizations of the structure functions. The 
K factor can be broken down to several characteristic components. K^i^-t accounts for the 
regularized virtual corrections, corresponding to the coefficient C; Kab B = g,q, q] for 
the real corrections as defined in eqs.(^OD. These K factors are shown for LHC energies 
in Fig. 11 as a function of the Higgs boson mass. For both the renormalization and the 
factorization scales, /i = M = rriH has been chosen. Apparently K^h and Kgg are of the 
same size and of order 50% while Kgg and Kgg are quite small. [Note that {K^rt + ^Kab) 
differs from [Ktot — 1) since the cross sections cxo are evaluated with different NLO and 
LO as values in the numerator and denominator.] Apart from the threshold region for 
Higgs decays into ti pairs, Ktot is insensitive to the Higgs mass. 

The absolute magnitude of the correction is positive and large, increasing the cross 
section for Higgs production at the LHC significantly by a factor of about 1.5 to 1.7. 
Comparing the exact numerical results with the analytic expressions in the heavy-quark 
Hmit, it turns out that these asymptotic solutions provide an excellent approximation 
even for Higgs masses above the top-decay threshold. For Higgs masses below ~ 700 
GeV, the deviations of the QCD corrections from the asymptotic approximation are less 
than 10%. 

There are two sources of uncertainties in the theoretical prediction of the Higgs cross 
section, the variation of the cross section with different parametrizations of the parton 
densities and the unknown next-to-next-to-leading corrections. Since all mass scales, 
the Higgs mass as well as the loop-quark mass, are very large, the notorious uncertainties 
from higher-twist effects can safely be assumed absent. 



^The dependence of the cross section on the factorization scale is very smah. 
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One of the main uncertainties in the prediction of the Higgs production cross section 
is due to the gluon density. This distribution can only indirectly be extracted through 
order Us effects from deep-inelastic lepton-nucleon scattering, or through complicated 
analyses of final states in lepton-nucleon and hadron-hadron scattering. Adopting a set 
of representative parton distributions |]63, |6^, B3] which are up-to-date fits to all available 



experimental data, we find a variation of about 7% between the maximum and minimum 
values of the cross section for Higgs masses above ~ 100 GeV, Fig. 12a. This uncertainty 
will be reduced in the near future when the deep-inelastic electron/positron-nucleon 
scattering experiments at HERA will have reached the anticipated level of accuracy. 

The [unphysical] variation of the cross section with the renormalization and factor- 
ization scales is reduced by including the next-to-leading order corrections. This is 
demonstrated in Fig. 13 for two typical values of the Higgs mass, mu = 150 GeV and 
ttlh = 500 GeV. The renormalization/factorization scale fi = M is varied as /i = prriH 
for p between 1/2 and 2. The ratio of the cross sections is reduced from 1.62 in leading 
order to 1.32 in next-to-leading order for itlh = 500 GeV. While for small Higgs masses 
the variation with for p ~ 1 is already small at the LO level, the improvement by the 
NLO corrections is significant at the NLO level for large Higgs masses. However, the fig- 
ures indicate that further improvements are required since the p dependence of the cross 
section is still monotonic in the parameter range set by the scale of order ttih- These 
uncertainties associated with higher-order corrections appear to be less than about 15% 
however. 

If the total energy is reduced from ^/s = 14 TeV to 10 TeV the production cross section 
for the SA4 Higgs boson decreases by a little less than a factor 2 for small Higgs masses 
and a little more than 2 for large Higgs masses. Fig. 12b. The K factors agree within less 
than ~ 5% for the two energies. 



3 The Neutral SUSy Higgs Particles 
3.1 The Basic Set-Up 

Supersymmetric theories are very attractive extensions of the Standard Model. At low 
energies they provide a theoretical framework in which the hierarchy problem in the Higgs 
sector is solved while retaining Higgs bosons with moderate masses as elementary particles 
in the context of the high mass scales demanded by grand unification. The minimal 
supersymmetric extension of the Standard Model {M.SSM.) |^ may serve as a useful 



guideline in this domain [70]. This point is underlined by the fact that the model led to a 



prediction of the electroweak mixing angle pS] that is in striking agreement with present 



high-precision measurements of sin 9w [29|. Although some of the phenomena will be 



specific to this minimal version, the general pattern will nevertheless be characteristic to 



more general extensions [p3, 71] so that the analyses can be considered as representative 



for a wide class of SUSy models. 
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Supersymmetry requires the existence of at least two isodoublet scalar fields $i and 



$2; thus extending the physical spectrum of scalar particles to five The AiSSJ^ 

is restricted to this minimal extension. The field $2 [with vacuum expectation value f 2] 
couples only to up-type quarks while $1 [with vacuum expectation value Vi] couples to 
down-type quarks and charged leptons. The physical Higgs bosons introduced by this 
extension are of the following type: two CP-even neutral bosons and [where 
will be the hghtest particle] , a CP-odd neutral boson A'^ [usually called pseudoscalar] and 
two charged Higgs bosons H"^. 

Besides the four masses m^o , rriHO, m^o and rriHi , two additional parameters define the 
properties of the scalar particles and their interactions with gauge bosons and fermions: 
the mixing angle (3, related to the ratio of the two vacuum expectation values tgjS = V2/V1, 
and the mixing angle a in the neutral CP-even sector. Supersymmetry gives rise to several 
relations among these parameters and, in fact, only two of them are independent. These 
relations impose a strong hierarchical structure on the mass spectrum [m/jo < mz,my^o < 



mijo and my/ < mH±\ which however is broken by radiative corrections [^, |3T| due to 
the large top quark mass. The parameter tg/3 will in general be assumed in the range 
1 < tg/9 < rrit/mh [7r/4 < /3 < 7r/2], consistent with the restrictions that follow from 
interpreting the M.SSM. as the low energy limit of a supergravity model. 

The M.SSM. Higgs sector is generally parameterized by the mass m^o of the pseu- 
doscalar Higgs boson and tg/3. Once these two parameters [as well as the top quark mass 
and the associated squark masses which enter through radiative corrections] are specified, 
all other masses and the mixing angle a can be predicted. To discuss the radiative correc- 
tions we shall neglect, for the sake of simplicity, non-leading effects due to non-zero values 
of the supersymmetric Higgs mass parameter and of the parameters At and Ab in the 
soft symmetry breaking interaction. The radiative corrections are then determined by the 
parameter e which grows as the fourth power of the top quark mass irtt and logarithmically 
with the squark mass M5, 

^ _ 3a 1 1 Yo (l iA<o) 

^ 27r s^c^ sin^ f5m% ml) 

with = 1 ~ = sin^ 6w [The main part of the two-loop effects can be incorporated 



by using the running MS top mass evaluated at the pole mass [[72| .] 

These corrections are positive and they shift the mass of the light neutral Higgs boson 
upward with increasing top mass. The variation of the upper limit on iriho with the 
top quark mass is shown in Fig. 14a for Ms = 1 TeV and two representative values of 
tg(3 = 1.5 and 30. While the dashed lines correspond to the leading radiative corrections 
in eq.®, 

mlo < m| cos^ 2/3 + e sin^ (3 (47) 

the solid hues correspond to the Higgs mass parameter fi = —200, 0, +200 GeV and the 
Yukawa parameters At = Af, = 1 TeV. The upper bound on m/iO is shifted from the 
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tree-level value mz up to ~ 140 GeV for rrit = 174 GeV. 

Taking m^o and tg/3 as the base parameters, the mass of the lightest scalar state 
is given to leading order by 



1 



m^o + rriz + e 

' (48) 



AO 



+ m| + e)2 — Am\om\ cos^ 2/3 — Ae{m\o sin^ /3 + m| cos^ /3) 



The masses of the heavy neutral and charged Higgs bosons follow from the sum rules 

2 2,2 2 I 

mjjo = rrij^o + rriz — rrif^o + e 

mjj± = m\o + m^y (49) 

In the subsequent discussion we will assume for definiteness that rrit = 174 GeV, Ms = 
1 TeV and = At = Aj, = 0. For the two representative values of tg/? introduced above, 
the masses 171^0,171^0 and mH± are displayed in Figs. 14b-d as a function of m^o. [The 
dependence of the masses on the parameters fi, At, Af, is weak and the mass shifts are 
limited by a few GeV 0.] 



The mixing parameter a is determined by tg(3 and the Higgs mass m^o, 

tg2a = tg2/j , r»^o +m| ^.^^ _ ^ < (50) 

mj^o — rn^ + e/ cos 2p 2 



The couplings of the various neutral Higgs bosons to fermions and gauge bosons depend 
on the angles a and (3. Normalized to the SAi Higgs couplings, they are summarized in 
Table |l]. The pseudoscalar particle A^ has no tree level couplings to gauge bosons, and 
its couplings to down (up)-type fermions are (inversely) proportional to tg/5. 







9<S>dd 


g^vv 


SM 


H 


1 


1 


1 


MSSM 




cos a/ sin /5 


— sin a/ cos (3 


sin(/3 — a) 






sin a / sin (3 


cos a / cos P 


cos(/3 — a) 




A' 


l/tg/3 








Table 1: Higgs couplings in the M.SS^A to fermions and gauge bosons relative to SAi 
couplings. 



Typical numerical values of these couplings are shown in Fig. 15 as a function of m^o 
and for two values of tg(3. The dependence on the parameters and At, Ab is very weak 
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and the leading radiative corrections provide an excellent approximation [^|. There is in 
general a strong dependence on the input parameters tg/3 and m^o. The couplings to down 
(up)-type fermions are enhanced (suppressed) compared to the SJ^ Higgs couplings. If 
m^o is large, the couplings of to fermions and gauge bosons are SAi like. It is therefore 
very difficult to distinguish the Higgs sector of the A4SSAi from that of the SMI, if all 
Higgs bosons, except the lightest neutral Higgs boson, are very heavy. 

Apart from cascade decays in some corners of the SlASy parameter space, the main 
decay modes of the neutral Higgs particles are in general bb decays [~ 90%] and r+r" 
decays [~ 10%], and top decays above threshold. The branching ratios for all the dominant 
decay modes are shown in Fig. 16. The gold-plated ZZ decays of the S^A Higgs particle 
above 140 GeV play only a minor role in the SUSy Higgs sector — and in large parts of 
the parameter space their role is even negligible. The total widths of the states remain 
small, 0{1 GeV), anywhere in the intermediate mass range and they do not exceed a few 
tens of GeV even for Higgs masses of the order of 1 TeV, Fig. 17. 

In addition to the conventional decays into SM. particles, the Higgs particles may also 
decay into chargino and neutralino pairs ||33|, Depending on the details of the SUSy 



parameters, the branching ratios for decays into these channels can add up to a few tens 
of percent; invisible LSP (lightest neutralino) decays, in particular, can even dominate in 
some domains of the M.SSM. parameter space. When kinematically allowed, the Higgs 
particles also decay into squarks and sleptons, with generally small branching ratios, 
though. For the present experimental bounds on non-colored and colored supersymmetric 
particles, see Refs. and [Q, respectively. 

The neutral Higgs particles will be searched for mainly in the decay channels r+r^ and 
77 at the LHC |ll8|, [1^. Large QCD backgrounds render the analysis of the dominating 
bb final states very difficult. Nonetheless, detailed feasibility studies have demonstrated 
that the hb decay channel may be accessible in associated Wh^ and ti/bbh^ events if 



a set of strong detector requirements is met 25 . 



3.2 The Two-Photon Decay Widths 

Similarly to the Standard Model Higgs boson, the precise prediction of the 77 widths of 
the SUSy Higgs particles is motivated by several points. This rare decay mode provides 
the most important signature for the search of the light Higgs bosons at hadron colliders. 
The values of the coupling constants are affected by the charged particle loops of the entire 
SUSy spectrum with masses far exceeding the light Higgs mass. The effect however is 
small in general for heavy SUSy particles since the main component of their masses is 
not generated by the Higgs mechanism so that these particles decouple asymptotically. 

The 77 coupling to Higgs bosons in supersymmetric theories is mediated by charged 
heavy particle loops built up by W bosons, standard fermions /, charged Higgs bosons 
H^, charginos c and sfermions / in the scalar cases h^,H^, and standard fermions and 
charginos [in the absence of sfermion mixing] in the pseudoscalar case A^. Denoting the 
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amplitudes hj Aj etc., the 77 decay rates are givenQ by 



r{n' ^ 77) 



128V27r3 



/ 



and 



r(A° ^ 77) 



32^2713 



j:Nce}9fAf + j:gtA 



(51) 



(52) 



The spin 1, spin 1/2 and spin amphtudes read to lowest order for the scalar Higgs bosons 



A^ 

^1/2 
A? 



_[2r2 + 3r + 3(2r-l)/(r)]/r2 
2[r + (r-l)/(r)]/r2 



- r 



and for the pseudoscalar Higgs boson 



At 



/2 



/(r)/r 



(53) 



(54) 



As usual, the scaling variable is defined as r = m|/4m^ with denoting the loop 
mass. The universal scaling function /(r) is the same as in eq.(|^). The coefficients gf 
denote the couplings of the Higgs bosons to W bosons, top and bottom quarks given in 
Table |I] and the couplings to sfermions and charginos which are recollected for the sake 
of convenience in Table ^ in the absence of sfermion mixing. [Including mixing effects in 
the scalar squark sector due to the soft parameters At, Ai, and /i does not change the 
production cross sections and photonic decay widths of the SUSy Higgs bosons in most 
of the parameter space, except in small regions where they play a significant role and lead 
to an enhancement of the signal |76[| .] 

Since the contributions of the squark loops are strongly suppressed compared to t, h 
loops, we shall restrict the discussion of the QCD corrections to the standard quark loops. 
These corrections will be parameterized again as 



A 



Q 



aLO 



l + C 



a. 



TC 



(55) 



The coefficient C depends on r = m'^ / ArriQ^fiQ) , where the running mass niQ^^Q) is 
defined at the renormalization point fiQ, 



C = Ci[mQ{fil)] + C2[mQ{fil)] log 



^4 



(56) 



^The scalar particles will generically be denoted by and all the neutral Higgs particles by 
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SM H 





MSSM 




MP - «) + 

cos(/3 -a)- 



2^(Sii cosa-Q,i sin«) 
2^{Sii sin« + g,, cos a) 
2^{-SucosP-QusmP) 





fL,R 


H 





MSSM 

A' 


^9'}T'-^{li - Cf Sin' e^)sm{a + (3) 

f s 

^fi ± ^(4^ - sin' ^1^) + /5) 




Table 2: Higgs couplings in the MSSM. to charged Higgs bosons, charginos and sfermions 
relative to SM. couplings. Qa and Sa (i = 1,2) are related to the mixing angles between the 
charginos ci and 62, Ref . ||70t| . 



The renormalization point is taken to be fig = m<j,/2; this value is related to the pole 
mass by the QCD formula noted in eq.(|^). The lowest order amplitude Aq^ must be 
evaluated for the same mass value mQ^fiq = [m$/2]^). The choice fig = m$/2 of the 
renormalization point ensures, a posteriori, a behavior of the 77 couplings which is well 
controlled for Higgs masses much larger than the quark mass. The QCD coupling as is 
evaluated at = m$ for A^.^ = 214 MeV. 

To regularize the pseudoscalar amplitude involving the 75 coupling, we have adopted 
the 't Hooft-Veltman prescription . A technical remark ought to be added on a subtle 
problem related to this implementation of 75 which reproduces the axial-vector anomaly to 
lowest order automatically. The multiplicative renormalization factor of the scalar {QQ) 
current is given by Z-yigg = 1 — ^2^^ where Z2, Zm are the wave-function and mass 
renormalization factors, respectively. To ensure the chiral-symmetry relation Y^ip'^-p) 
75r(p',p) in the limit mg for the fermionic matrix element of the pseudoscalar and 
scalar currents, the renormalization factor of the pseudoscalar current has to be chosen 
W\ as 



ZaQQ — Z-HQQ + n — (57) 
O 77 
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The additional term, supplementing the naive expectation, is due to spurious anomalous 
contributions that must be subtracted by hand. 



The Limit of Large Loop Masses 

For large rriQ, the coefficient C2 in eq. (|56D is of order l/rriQ and approaches zero for the 
scalar and pseudoscalar Higgs bosons. It has been shown before that ci approaches (— 1) 
for scalar Higgs bosons; for the pseudoscalar Higgs particle ci vanishes asymptotically, i.e. 

niQ ^ oo : ^ -1 

cf ^ (58) 

This result for the pseudoscalar Higgs boson can also be derived from the non-renormali- 
zation of the anomaly of the axial-vector current. In the same way in which the Ti.'-ff 
coupling in the local limit can be related to the anomaly of the trace of the energy- 
momentum tensor, we can derive the ^4^77 coupling from the anomaly of the axial-vector 
current ||79| , 

d^Jl = 2mQQi^5Q + N,e%—F^,F^, (59) 

with F^y = tfj^yapFafj denoting the dual field strength tensor. Since, as famihar from the 
Sutherland- Veltman paradox, the matrix element (77|9^j5 |0) of the divergence of the 
axial-vector current vanishes for zero photon energy, the matrix element (77|w^qQ^75<5|0) 
of the Higgs source can be linked directly to the anomalous term in eg . (|59|) . It is well- 
known that the anomaly is not renormalized if the QCD strong interactions are switched 
79[. As a result, the effective ^4^77 Lagrangian 



on 



C,ff{A^ll) = N,el^ (V2Gf)"^ F,J,,A'' (60) 



in 

is valid to all orders of perturbation theory in in the limit m\o <^ 

The Limit of Small Loop Masses 



Also in the opposite limit of small quark-loop masses compared with the Higgs masses, the 
'H°77 and y4''77 couplings can be calculated analytically. This limit is useful in practice 
for large tg/3 values where the b quark coupling to the heavy Higgs bosons and is 
strongly enhanced. As anticipated theoretically, the leading and subleading logarithmic 
terms are chirally invariant and we obtain the same QCD correction in this limit for the 
scalar and pseudoscalar couphngs, 

mgifiD^O : C^'^-.-llog2(-4r-ze)-|log(-4r-ze)+21og4- (61) 

The finite non-logarithmic contributions to C may be different in the scalar and pseu- 
doscalar cases. 

The amplitudes for scalar loops and for pseudoscalar loops are shown in Fig. 18 
as a function of r The coefficients are real below the quark threshold r < 1, and 
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complex above. Very close to the threshold, within a margin of a few GeV, the present 
perturbative analysis^ can not be applied. [It may account to some extent for resonance 
effects in a global way.] Since QQ pairs cannot form 0+"*" states at the threshold, QmC-j^ 
vanishes there; ^eC-n develops a maximum very close to the threshold. By contrast, since 
QQ pairs do form 0~"^ states, the imaginary part Q^mC^ develops a step that is built up 
by the Coulombic gluon exchange [familiar from the Sommerfeld singularity of the QCD 
correction to QQ production in e~^e~ annihilation]; ^cCa is singular at the threshold. 

The singular behavior of the A^'j'-f coupling demands a more careful analysis at the 
quark threshold The form factor is given to lowest order near the threshold by 



K^'^'irQ) = f{rQ)/rQ - ^ + for tq ^ 1 (62) 



Where (3 = Jl — Tq^ is the quark velocity above the threshold. The QCD corrections to 



the imaginary part can be found by attaching the Sommerfeld rescattering correction [80| 



Ccoui = — ^ ^l + -Z for Z = ^ (63) 

which corresponds to the exchange of a ladder of Coulombic gluon quanta between the 
slowly moving quarks. The QCD corrected imaginary part of the ^4*^77 coupling may thus 
be written 

2 

'^m = irpCcoui = 7^13+ -n^Os (64) 

approaching a non-zero value at threshold. The real part can be derived from a once- 
subtracted dispersion relation so that near the threshold 

4 ^ ^'^'^ + ^ [- log(rQ - 1) + ^vr + const] (65) 

The smooth constant term needs not be fixed if we analyze only the singular behavior. 
For the QCD correction near the threshold we therefore obtain the simple relations, 

Tq ^ 1 : ^eC^ ^ \og{TQ - 1) + const 

g 

^ +-TT ^ 8.38 (66) 
3 

The absolute size of the imaginary part and the logarithmic singularity of the real part 
are in agreement with the numerical analysis presented in Fig. 18b. 

In Figs. 19a,b the QCD corrected 77 widths for the h^, H^, Higgs bosons are dis- 
played, taking into account only quark and W boson loops for two values tg/3 = 1.5 and 



^By choosing the renormahzation point jiQ — m<s>/2 the perturbative threshold Eth = '^nT-Qi^g) 
coincides with the on-mass shell value proper. A shift between m$/2 and m$, for instance, affects the 
widths very little away from the threshold. 
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tg/5 = 30. While in the first case top loops give a significant contribution, bottom loops 
are the dominant component for large tg(3. The overall QCD corrections are shown in 
Figs. 19c, d. The corrections to the widths are small, ~ 0{as/7^) everywhere. [Artificially 
large values of 6 occur only for specific large Higgs masses when the lowest order ampli- 
tudes vanish accidentally as a consequence of the destructive interference between W and 
quark-loop amplitudes, see also [Q.] Thus, the QCD corrections are well under control 
across the physically interesting mass range if the running of the quark masses is properly 
taken into account. 



3.3 The Gluonic Decay Widths 

The gluonic decays of the Higgs bosons 
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are mediated by quark and squark triangle loops. In the same notation as in the preceding 
section we find for the widths in lowest order 

2 



r{h' - 99) 



r(i/° 99) 



r(A° gg) 



7T^ 



7 Yl 9qAq + jY1 9^A^q 
^ Q ^ - 



Gpa 



2 



Q 



Q Q 



1672' 



-m" 



E9^A^ 



Q 



(67) 



Since the contribution of heavy squark loops is small, we will neglect these effects in the 
following discussion and we will focus on the dominant quark contributions. 

The QCD corrections to the gluonic decay widths are large. Besides the virtual correc- 
tions, the widths are affected by three-gluon and gluon plus light quark-antiquark final 
states, 

h'^^H'^A'-^ggg and gqq (68) 

Proceeding in the same way as for the Standard Model, the result can be written in the 
form [$ = /lO, i/o, AO] 



r($^ (?(?((?), 



^Loi"^ 99) 



i+^*(t-q: 



a. 



vr 



(69) 



with 



^ / X 95 7,, 33-2Nf , 

4 6 6 mi 



Ea{tq) 



''H 
,2 



97 7 33 - 2Np ^ fi- 

- -Nf + — - log -V + ^Ea 



4 6 



6 



(70) 
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In the limit of large loop masses, a contribution 11/2 to the coefficients for scalar 
states is related to the effective Lagrangian after the heavy quarks are integrated out; the 
remaining part is associated with the rescattering and splitting corrections. As a result 
of the non-renormalization of the axial anomaly, the coefficient for the pseudoscalar state 
is entirely due to the rescattering and splitting corrections. The corrections A£'$ are dis- 
played in Figs. 20a,b as functions of the corresponding Higgs masses within their relevant 
mass ranges for tg/3 = 1.5 and 30. Due to the bottom contribution, the deviations from 
the heavy quark-loop limit are significantly larger than in the SAi case, thus rendering 
this limit useful only for tg(3 close to unity. In Fig. 21 the gluonic decay widths including 
the QCD radiative corrections are presented for tg(3 = 1.5 and 30. They are enhanced by 
about 50% to 70% as a result of the large QCD corrections. In a margin of a few GeV near 
the threshold [mA ~ 2mf] the perturbative result of the pseudoscalar decay width is not 
valid due to the Coulomb singularity in analogy to the photonic decay 77. The 

final branching ratios of all decay processes in the AiSSAi are shown in Fig. 16. For the 
light Higgs particle the gluonic decay mode is significant only for masses close to 
the maximal value, where has SAi like couplings. For the gluon decay mode is 
significant only slightly below the top-antitop threshold and for small values of tg/3 where 
the coupling to top quarks is sufficiently large. For the pseudoscalar Higgs boson A^, 
the gluonic decay mode is important for small values of tg/5 and below the top-antitop 
threshold, where it can reach a branching fraction of ~ 20%. 

In the limit of large quark masses, the Higgs-gluon-gluon coupling can be described 
by gauge-invariant effective Lagrangians, 

r - ^ f for \ /^("«) r"- -wO 



I^A^,, = ^ {^Gf) G%G%A<^ (71) 



with f3 and 7^ defined previously. They take account of the local interactions of the 
particles involved and serve as kernels for the standard gluon and light-quark corrections. 

3.4 Higgs Production in pp Collisions 

The production of SUSy Higgs particles at hadron colliders has received much attention 



in recent years after the pioneering investigations in Refs. |3^. The situation is critical 
since the ffist analyses could not ensure that the entire AiSSAA. Higgs parameter space 
could be covered at the LHC. Yet, high statistics analyses appear to solve this problem 



if the decays to SAi particles are dominant |T^. A second similarly severe problem 



has arisen from the difficulty to detect the heavy Higgs particles for masses above a few 
hundred GeV and moderate values of tg/3 where the production rates in the experimentally 
clear t^t~ channel are too small to be exploited in practice. However, no final picture 
has emerged yet, since the detailed conclusions depend strongly on the detector design. 
Additional decay and production channels, based on the tagging of heavy quarks, may 
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also help close the hole in the parameter space |7H . 



The dominant production process for SUSy Higgs particles at the LHC is the gluon 
fusion mechanism. Besides the virtual corrections, the bremsstrahlung of additional glu- 
ons, the inelastic Compton process and quark-antiquark annihilation, 
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h^lH^A^g) and gq^ hynyA^q, qq^h^'/HyA'g 



/ uO / aO, 



,0/rjO /4O, 



contribute to the Higgs production. The diagrams relevant to the various subprocesses 
are the same as for the Standard Model in Fig. lb. The parton cross sections may thus 
be written 



1 + C{tq 



a. 



TT 



<5(l-f) + A,(r,rQ)^e(l-f)| 

TT J 



(72) 



for i,j = g,q,q and f = m|/s. The final result for the pp cross sections can be cast into 
the compact form [$ = H^, A^] 



a{pp ^ h^/H^/A^ + X) = do 



C 



TT 



r<s>— h Aagg + Aagg + AcTqq 



(73) 



with 



h,H 
(^0 



288^271 
12872' 



4 H 9qA^ 



TT 



T.9$A^ 



(74) 



after folding the parton cross sections with the MS renormalized quark and gluon densities 
[tq = m'^/ArriQ and r<j, = m^/s]. The virtual/IR and hard corrections have the same 
generic form as before, eqs. (|5P|-^TD. As a result of the factorization theorem, the parity 
and the specific couplings of the Higgs bosons are not relevant for the infrared/coUinear 
form of the cross sections, related to interactions at large distances. The specific properties 
of the Higgs bosons affect only the non-singular coefficients c and d in eqs. (p9|-p|) . 

In the limit of large quark-loop masses compared with the Higgs masses, only the 
coefficients c depend on the parity of the Higgs particle. 



Tq = ml /Ami 







11 

y 

6 



(75) 



The coefficients d are universal. The next-to-leading term in the expansion for the scalar 
Higgs bosons has also been calculated analytically ||^ . The form of the cross sections for 
the parton subprocesses in the heavy quark-loop limit if the final states are analyzed, is 
given by the same expressions as eq. (43). In the opposite limit of small quark-loop masses. 
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chiral symmetry is restored for the leading and subleading logarithmic contributions to 
the coefficients, which are given by the same expressions as eq. (44). 

The final results of the pp cross sections are predicted in the subsequent figures for 
the LHC energy ^/s = 14 TeV. [A brief summary is appended for = 10 TeV.] Again, 
the two representative values tgf3 = 1.5 and 30 arc chosen and the top mass is fixed to 
rrit = 174 GeV. If not stated otherwise, we have adopted the GRV parameterizations of 
the quark and gluon densities. For the QCD couphng we have chosen the average value 
a^^\mz) = 0.117 in the final cross sections while the discussion of the K factors is carried 
out consistently in the GRV frame. [The GRV NLO coupling is close to the lower la 
boundary of the global as fit.] 

The K factors, Ktot = o'ho/'^lo: ^-re defined by the ratios of the HO cross sections 
to the LO cross sections. They are shown for LHC energies in Fig. 22. They vary little 
with the masses of the scalar and pseudoscalar Higgs bosons in general, yet they depend 
strongly on tg/3 as shown in Fig. 23. For small tg/3, their size is about the same as in 
the SA4, varying between 1.5 and 1.7; for large tg/3 however they are in general close 
to unity, except when approaches the SAi domain. The cross sections are shown in 
Fig. 24. Apart from exceptional cases, they vary in the range between 100 and 10 pb for 
Higgs masses up to several hundred G^. Beyond ~ 300 they drop quickly to a level 
below 10~^ pb. Similarly to the SM., a factor of about 2 is lost if the pp coUider energies 
is reduced to 10 TeV, Fig. 25. 

The variation of the cross sections with the renormalization/factorization scale is re- 
duced by including the next-to-leading order corrections. The dependence of the cross 
sections for low masses. Fig. 26, is of order 15%; the fi dependence remains monotonic. 
Thus the next-to-leading order corrections stabilize the theoretical predictions for the 
Higgs particles in the intermediate to large mass range, yet further improvements must 
be envisaged in the future. 

It is apparent from the previous figures that the next-to-leading order corrections 
increase the production cross sections for the SUSy Higgs particles, in some areas of the 
parameter space even strongly. 

4 Summary 

We have presented a complete next-to-leading order calculation for the production of 
Higgs particles at the LHC in the Standard Model of the electroweak interactions as well 
as in its minimal supersymmetric extension. These corrections stabilize the theoretical 
predictions compared with the (ill-defined) leading-order predictions. The QCD radiative 
increase the production cross sections significantly so that experimental opportunities to 
discover and detect these fundamental particles increase. 
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APPENDIX A: The 7^77 and ^77 Couplings 

In this Appendix, we summarize the complete analytical result for the QCD corrected 
CP-even 7-^77 and CP-odd ^77 vertex form-factors, in the case of arbitrary Higgs boson 
and quark masses. 

As discussed in sections 2.1 and 3.1, the radiative QCD corrections to the quark 
contribution to the two-photon Higgs boson decay amplitudes can be written as 



_ aLO 

where the coefficients C$ split into two parts. 



1 + C^^ 



(A.l) 



m 



(A.2) 



Q 



with the two functions Cf and depending only on the scaling variable tQ, 

r = m|/(4mj) = p$/4 (A.3) 
For the CP-even and CP-odd Higgs bosons, the coefficients Cf and are given by 

C? = 



- 2r-iFo^(r)] " {2t~\8t~' - 17) + 4t-\3t-' - 4)/(r)(7(r) 
+4r"2(^-2 _ 2^-1 + 3)j(^) + 24r"2(_^-i ^ 2)^(r) + 2r"^(-5r-^ + 6)/(r) 



+r-'{-2r-' + + 3)A;(r) + Qt-^t'^ - l)h{r) - 8( 
+Ar-\2r-' - 3) [/a -h + h] + 16r-\r-' - 1)1, } 



_-2 



3r 



-1 



2)/i 



-3r 



2r-i-l)/(r)-(7(r)-l} 

3r-iFo^(r)] ' {24r-' + 16r-V(r)^?(r) + 6T-\r-' - 2)/(r) - 36r-2^(r) 
-12r-2/(r) - 3r-2(r-i + l)A;(r) + 6T~^h{T) + 8(-r-i + 2)/i 
+ 12r-i[/2-/3 + /4] + 16r-i/5 } 



with 



2r 



-1 



l + (l-r-i)/(r)l and ^^^(r) = r-i/W 



In terms of the auxiliary variables 



a± = (1 ± Vl - r-i)/2 



(A.4) 



(A.5) 



(A.6) 
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Singularities are fixed by attributing to the quark-loop mass a small imaginary part: itiq —f tjiq — it 
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the functions /, /, k and h are defined by 

r > 1 



1{t) 



arcsm \/t 
1 



log in 



Vr-i - 1 



arcsm v r 



— 1 



log 27r 



Li3(l/a+) +Li3(l/«_) 



r < 1 
T > 1 



fc(r) = 5i,2(l/«+) + 5i,2(l/«-) 



Mr) 



+ 2 



+ 2C(3) (A.7) 



Here, Li2, Li^ and 6*1,2 are polylogaritlims, defined |77 



as 



Lz2(x) = / — log(l-xy) 
Jo y 



Li3{x) = / — log?/log(l - xy) 
Jo y 

Si,2{x) = \f'^\oi{l-xy) 
I Jo y 



(A.8) 



The expressions of /i....,5, which have been reduced from four- and five-dimensional down 
to one-dimensional Feynman integrals, are much more involved: 



p I — px[l — x) y \l — px ^ 



px 



1 — px{l — x) 



-Li2 \^ 
]-Lt2 



1 — px j 
/ —px{l — x) 
— px{l — x) 



log[x{l — x)] log[l — px] 



1 - ^ 



+ 



log^(l — px) 



dx 



1 - 



+ log(l - px) log 



1 - ^ 



px 



1 — px{l — x) 



p[l — px{\ — x)\ 



Li3[px(l - X)\ + LZ3 I I - -^^^3 



-px 



-px" 
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+51,2(1 -x)- ^i,: 



1 — X 



-X 



(l-px)) 



/ 1 — px \ lo^^ X 
-Li2 { — 1 log a; H — log(px) 



+Li3 



-px 



-X 



1 — X 



-X 



(1 — x){l — px) 
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log(l - px)Li2 
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log^(l - px) 



(A.9) 



1 dx 



px 6 — log ( 1 — ) log ( 1 — 
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1 — px{l — x) ) 
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In terms of the variables a = {ad — hc)/d and (5 = b/d, and the functions 
F,{a,b) 



1 /^l + «\r- fbl + a\ 1, 2/l + «\, / 
+ log {—^) L,, + - log^ (_j log ^ 
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the expressions K^, which appear in the integral /2, are given by [K^ and will be used 
later onl 
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1 + a 



1 + log(l + g) 
1 + a 



+ 



a 

a — b 



log(l + a) {log(l + 6)-l} 



ab 



log 



log(l + b) + Li2 



1 + b 
b 

—a 
) — a 



log(l + &) 



— Lio I —a 



1+6' 



Ki,{a,b) = ^{^log(l + 6)-i±^log(l + a)} 



(A.16) 



APPENDIX B: The Infrared Regularized Virtual Corrections for 
the Tigg and Agg Couplings 

The complete analytical expressions for the virtual, infrared regularized, QCD radiative 
corrections to the Tigg and Agg couplings are summarized in this appendix. As discussed 
in section 2.3 and 3.3, the virtual corrections split into an infrared part tt^, a logarithmic 
part depending on the renormalization scale n and a finite piece depending on r 



C-7r2 + log^ 

D m| 

Again the coefficient c$ can be split into two parts 



+ c$(r) 



(B.l) 



(B.2) 



where is the scale at which the quark mass is defined; the coefficients 5i, B2 read in 
the case of CP-even and CT'-odd Higgs bosons, 

-1 



Bf 



16r-^Fo^(r) 



H 



B. 



Bt 



{32r-2(-3r-i + l)/(r)^(r) 

+40t-2(-2t-^ + 4r-^ + 3)/(r) + 192r-2(4r-^ - 5)^(t) 

+216r"^(-2r"^ + 3)r(r) + 36r~2(-2r-^ + 3)p(r) 

+4r-^(18T-^ - 49r-^ + 30)/(r) + 2t-^(38t-^ - 199t-^ + 186)/c(r) 

+3t-^(-18t-^ + 65t-^ - 53)/i(r) + 144(r-^ - l)[2/6 - r/g] 

+16t-^(-29t-^ + 20) + 128(-T-^ + Sr"^ - 2)h 

-8(2t-^ - 3)[t-^72 - lOr-^/s + r'^h + I8/7] - 32t-^(t-^ - l)h } 

-6r-i [Fo^(r)]"' {(2r-^ - l)/(r) - ^(r) - l} 

'2At-^F^{t)[' {-32r-V(r)^?(r) - 120r-2(r-i + l)/(r) + 1152r-2^?(r) 
-648r-V(r) - lOSr^Vl^) + 12r-2(9r-i - 10)/(r) + 6r"2(19r-i - 62)A;(r) 
+3r-2(53 - 27r-^)/i(r) + 128(2 - r-^)/i - 24[t-^/2 - lOT-^/3 + r'^h + I8/7] 
-32r-i/5 + 144[2/6 - r'^h] - 696t-2 } 



Bi, 



4t-i 



F,\t)\-" {/(r)-^(r)/(r-^-l)} 
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(B.3) 



The functions , are given in eq. (A. 5), while the functions /, Z, k and h are given 
in eq. (A. 7); the two remaining functions r and p read 



r(T) = --Vl -T- 



Li-2 



q;_|_ — q;_ 



— Lio 



+ ^log[4(l-r) 



p(r) = Vl - r-i [Li2(l/a+) - Li2(l/a_)] 



(B.4) 



The integrals /i to /5 were presented in eqs. (A.9-A.13); the integrals Je, h and /s are 
given by 

9 /"l rlr (R 

h = - 1 . {|V[l-p^(l-^)]-21og(l-x)log^[l-px(l-x)] 



p Jo — x) 

+31og[l - px{l - x)]Li2{px{l - x)) + 8^i,2(pa:(l - x)) 



+2 



px 



1 — px 



Ki —px{l — x), —X, 



px 



1 — px 



- Kn 1, -px{l - x), 1, 



px 



1 — px 



+4pa;(l — x)K3 1, —px(l — x),l — px, px^j | 



(B.5) 



2 dx ( 8 

^7 = - / — ^log'[l-px(l-x)]-21og(l-x)log2[l-px(l-x)] 
p Jo p 16 

+31og[l - px(l - x)]Li2(px(l - x)) + 85i,2(px(l - x)) 



+2 



px 



1 — px 



K4 —px{l — x), —X, 



px 



1 — px 



- Kn 1, -px{l - x), 1, 



px 



1 — px 



+4px(l — x)K3 1, —px(l — x),l — px, px^j I 

+ -J^dx xil - X) 1;^^^ [log[l - Px{l - x)]U2 ^— 

—px{\ — x) 



—px{l — x) 
px{l — x) 



-Si: 



1 — px{l — x) ^ 

+F2 [1, -X, 1, -px{l - x)] + F2 1, -px{l - x),l - px, px^ 

—F2 1, —X, 1 — px, px^ + Ki [1, —X, 1, —px{l — x)] 

—Ki 1, —px{l — x),l — px, px"^ — Ki 1, —X, 1 — px, px"^ 
-2K3[1, -px{l - x), 1, -x] - [1, -X, 1, -px{l - x)] 



-\-Kyj 1, —X, 1 — px, px' 



-K 
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1, —px{l — x),l — px, px^ 



-K18 [1, -px{l - x)] + K18 1 - px, px^ } 



+ - / da; a; < logfl — pxil — x) 
P Jo I 



px 1 — px 
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+ log^[l — px{l — x)] 



3 1 3 

+ -7Z T + 



2px p{l — x) 2(1 — px) 



+ log 



1 — px{l — x) 
I- px 



^2 — X ^1 — X 
-2 ^ + 2 



px^ 



X 



^1— Xl — px^ 2 ~ ~ ^) 



2x px 



1 — px 



px 



+ 2^ \ogx 



1 — px{l — x) ) 1 — px 



+ 



1 px 



■ log^ X — F2 1, —px{l — x),l — px, px^ 



+ F2 



1, —X, 1 — px, px^ 



21- px 

+Ki[l, —X, 1 — px, px"^] + Ki[l, —px{l — x),l — px,px'^] 

+2(1 - x)K2[l, -px{l -x),l- px,px^] + 2K2[l, -X, 1, -px{l - x)] 

+2pa;(l — x)K2\i, —x, 1 — px, px"^] 

p.r{l — x) 



-2- 



[K2[l, -(1 - x), 1, -px{l - x)] + K2[l, -px{l - x), 1, -(1 - x)]] 



1 — px 

+2(1 - x)K2[l, -px{l - x), 1, -x 



-2 



1 — px{l — x) 

(1 — pxy 

1 — px{l — x] 



Ka -X, 



px 



1 — px 



+ Kq -px{l - x). 



px 



1 — px 



(1 - pxy 

1 — px{l — x) 

(1 — pxY 



Ku I -X. 



px 



Kl2 -X, 



1 ~ px ^ 
px \ 



Kii -px{l - x), 



K12 -px{l - x) 



px 



1 — px 



px 



1 — px 



\ — px ^ 

—Kyj\\-, —X, 1 — px, px^\ + Kij[l, — pa;(l — x),l — px, px^] 
—Ki^\\ — px, px'^] + [1 — px{l — a;)]A'i9[l — px, px^] 

-[1 - px{l - x)] K2q[1, -px{l -x),l- px, px^] - K2q[1, -X, 1 - px, px^jj I (B.6) 



2 fir f 1 

h = - —\--\og{l-x)\og\l-px{l-x)]+2>S,,2{px{l-x)) 
p Jo X l 2 

-px{l — x) \ 



2 r^dx ( 1 
P 

—2 log[l — px{l — x)]Li2 



1 — px{l — x) 



—px{l — x), -{1 — px) 



X 



-px{l-x). 



1 ^ X 



X 



+px{l - x)K^ [I, -px{l - x),x, (1 - a;)(l - px)]} (B.7) 

While the functions Ki-Ku arc defined in eqs.(A.16), the remaining functions follow from 

1 r 1 + a 



= ^{YTilog'(l + «)-2«^2(l,a,l,6)| 



Ki^{a,h,c) = 



log(l + a) log(l + h) I \ a 



+ - 



a — c 



c(l + c) c ya — c 

-K2{l,h,l,c) + -—K2{l,a,l,h) 



K2{l,h,l,a) 



ac ,^ , _ . 6^ , be 



b-c 
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b-c 



K2{l,a, l,c) 



KM 
KiG{a,b) 
Ki7{a, b, c, d) 

KiQ{a,b) 
^20 (a, b, c, d) 



a - (1 + a) log(l + a) 

-(1 + a) log^(l + a) - 2aLi2(-a) 

Kr{a + b, -b) 

Ki{a + b, —b, c + d, —d) 



M- -C(2) 



\ a . 



2Lis I — 



ab ^ i a + 6 , 



b{a + b) 



Lio — - 



C(2) 
a(a + b) 



\\b_ 
d\ad 
1 



Lin 



-M 



a-\-b ^ 
b \ 



— Lio 



d 



c + d , 



c + d \a + b I c + d 



K2{a, b, c, d) 



(B.8) 



APPENDIX C: The Real Corrections for pp ^ A and 7Y, A ^ gg 

Finally, we give here the complete analytical expressions for the real corrections to the 
processes pp — > T-i/A and H/A ^ gg. We start with the corrections to the production 
process and define the variables 



T 



t 



m 



Q 



m: 



Q 



m: 



2 ' 



U 



u 



(C.l) 



m 



H 



i = — <§(! — f)^ , 



u 



-s{l - f){l - v) 



(C.2) 



n = -r 



T 
4 ' 



U 



T,, = — 



(C.3) 



The coefficients dqq,dgg and dgg which appear in the real QCD corrections, eq. (40), for 
the Higgs production, can be cast into the form 

2 



dgq{r,TQ) 



2.P 2 



EQi^O*(TQ) 

3 3 Jo V ^ 

l + (l-f)2(l-T;)2 



Q 



1 - 2 



1-f 



2^ 



dgg{r,TQ) 



1 — f Jo V 



1 dv 
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- V (T) 

-l-f^-(l-f)^ 



(C.4) 



with Fq given in eq. (A. 5). With the functions / and g given in eq. (A. 7), and A^^^ 
can be represented as 

1 fi 

•A%MT,U) = \Cf\'+\C!\' + \C!\' + \Ct\' (C.5) 
The functions Cf follow from 

CfiS, T, U) = -CfiT, S, U) Cl{S, T, U) = CfiU, T, S) (C.6) 
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(C.7) 



The coefficients T,- read 



Ti = 1 , T2 = 2f{TQ) , r3 = 2/(r,) , T, = 2f{n) {C.i 
Ts = 2/(r„) , r6 = 2[l-^(rQ)] , T, = 2[l - g{r,)] , Tg = 2 [1 - (7 (n)] 
Tg = 2[l-g (r„)] , Tio = J{S, T, U) , Tn = J{S, U, T) , T12 = J(T, 5, f/) 



with 



J(5,T,f/) = /3(5,T,f/,5) + /3(5,T,f/,f/)-/3(5,T,f/,p) 

/5+ 



h{S,T,U,X) 



1 2 L. ^ /3 

LZ2 



+ Li 



SU (3+ - (3. 

/9- 



/5_ — 



Lin 



3+ - a+ 



(C.9) 



/?_ — a+ 



— Li-; 



/3, 



[3- 



/9+ — a_ 



and 



and l±Jl + |^ 



(C.IO) 



The coefficients Pij for CP-even neutral Higgs bosons are given |Q by 

UT-S^ 



Pi, 



-12 ^ 



(f/ + 5)(T + 5) 
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A,3 



- -3{S- 4) 

4U^T + W^TS - U^S^ + AU^S + AUTS^ + 8US^ + 3"^- AS^ 



-3 



1 o TrT 



s{T + sy 

U^T + 2U^S + UT^ + 4UTS + 5US^ + 2T^S + 5TS^ + 4S^ 



^1,6 - - 



Plj 


= 


Pl,9 


= 12UT {T + 2S){T + S) 


PlM 


= 3TS{4-S)/2 


Paj 


= 12 p 


^4,3 


= 3(4-p) 


Pi,b 


= 3(4-p) 


Pi,7 


= 




= 


Pa,ii 


= 3TS{4-p)/2 



Similarly for CP-odd neutral Higgs bosons, 



Pl,l 


= 


Pl,2 


Pi,s 


= -2S 


PiA 


^1,6 


= 




Pl,7 


= 


Pi,s 


Pl,9 


= 


Pi, 10 


Pi,u 


= -TS^ 


Pi, 12 


P^l 


= 


Pa,2 


^,3 


= -2p 


Pa,a 


-^4,5 


= -2p 


Pa,6 


Pa,7 


= 


Pa,s 


Pa,9 


= 


Pa,io 


Pa,ii 


= -pTS 


Pa,12 



{u + sy{T + sy 

Pi,8 = 12 UT {U + 2S){U + 5) -2 

Pi,io = 3USiA-S)/2 

Pi 12 = -3 UTS-^ {AUT -S^ + 12S) /2 

Pa,2 = -9 (4 - p) 

P4,4 = 3(4-p) 

Pa,6 = 

^4,8 = 

P4,io = 3US{A-p)/2 

Pa,i2 = 3UT{4-p)/2 



= -2S{UT -US-TS- 3S^){U + S)-^{T + 5)-^ 

= 2S {U - S){U + S)-^ 

= 25(T-5)(r + 5)-^ 

= 

= -us^ 

= STU 

= 6p 

= -2p 

= 

= 

= -pUS 

= -pUT 
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For the QCD corrections to the gluonic decays of the Higgs bosons the correction 
factors A£'$ defined in eqs. (23) and (70) can be written as 



AE* , + AE* + NpAE'^ 



911 



11 



ca{tq) - 6 



A£; 



999 



AE. 



911 



I dxi [ 
Jo Ji- 



dX2 . ^ 
XI I StU 



mlA%g{S,T,U) mj + s^ + t^ + u^ 



2stu m| 



dx' 



s^ + v? \m% EQvA*g(T) 



2— 2~ 



(C.ll) 
(C.12) 
(C.13) 

(C.14) 
(C.15) 



EQi^o*(rQ) 

with ct-c and ca given in eqs. (B.2, B.3); the kinematical variables are defined as 

s = m^(l - X3) , i = m^(l - X2) , w = m^(l - xi) , xi + X2 + = 2 (C.16) 
and 



t 



,2 ' 



u 



rriQ rriQ ml 

for the Mandelstam variables normalized by the quark mass. 



(C.17) 
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Figure Captions 



Fig. 1: Generic diagrams of the gluon fusion mechanism gg ^ H for the production of 
Higgs bosons: lowest order amphtude (a), and QCD radiative corrections (b). 

Fig. 2: Generic diagrams for the amphtude of the Higgs boson decay into two photons 
H — > 77: (a) lowest order ly-boson amplitude, (b) lowest order quark amplitude, 
and (c) QCD radiative corrections to the quark amplitude. 

Fig. 3: Generic Feynman diagrams for the amplitude of the Higgs boson decay into gluons 
H —>■ gg: (a) lowest order amplitude, and (b) QCD radiative corrections. 

Fig. 4: The real and imaginary parts of the lowest order amplitudes Af (a) and Aw (b) of 
the i777 vertex as a function of Tf^w = 

Fig. 5: Real and imaginary parts of the radiative correction factor to the quark amplitudes 
for the if 77 coupling; the renormalization point for the quark mass is taken to be 
AiQ = ruH/^ in (a) and /xq = mg in (b). 

Fig. 6: (a) The QCD corrected partial decay width of the Higgs boson to two photons as a 
function of the Higgs mass, and (b) the size of the QCD radiative correction factor 
(in %). 

Fig. 7: (a) Comparison of the size of the infrared regularized virtual QCD corrections to the 
quark amplitude for the pole mass mQ{mQ) and the nmning quark mass mQ(m///2); 
for large quark-loop masses the coefficient C approaches the value C = 7r^/2 + ll/4; 
(b) The deviation /S.E of the radiative QCD correction to the decay H ^ gg from 
its value in the heavy quark-loop limit; the renormalization scale is taken to be 
H = ruH- 

Fig. 8: (a) The QCD corrected partial decay width of the Higgs boson into two gluons 
(in MeV) as a function of the Higgs mass, and (b) the size of the QCD radiative 
correction factor; the renormalization scale is taken to be = ruH- 

Fig. 9: (a) Total decay width (in GeV) of the Standard Model Higgs boson as a function 
of its mass, and (b) the branching ratios (in %) of the dominant decay modes 
{mt = 11 A GeV). All known QCD and leading electroweak radiative corrections are 
included. 



Fig. 10: (a) Feynman diagram for the effective couplings of the Higgs boson to gluons in the 
heavy-quark-loop limit, and (b) generic Feynman diagrams of the effective QCD 
corrections to the decay H ^ gg in the heavy-quark-loop limit. 

Fig. 11: K factors of the QCD corrected cross section a{pp H + X) aX the LHC with 
cm. energy ^/s = 14 TeV. Kmrt and Kab (A, B = q,g) are the regularized virtual 
correction and the real correction factors, respectively; Ktot is the ratio of the QCD 
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corrected total cross section to the lowest order cross section. The renormalization 
and factorization scales are taken to be /i = M = niH and the GRV parameteriza- 
tions for the parton densities have been used. 

Fig. 12: (a) The spread of the Higgs boson production cross section at the LHC with cm. 

energy of ^/s = 14 TeV for two parameterizations of the parton densities, (b) The 
total Higgs production cross section at the LHC for two different cm. energy values: 
^/s = U TeV and = 10 TeV. 

Fig. 13: The renormalization and factorization scale dependence of the Higgs production 
cross section at lowest and next-to-leading order; the Higgs mass is chosen to be 
(a) ruH = 150 GeV, and (b) tuh = 500 GeV. 

Fig. 14: (a) The upper limit of the lightest scalar Higgs boson mass in the AiSSM. as a 
function of the top quark mass for two values of tg/5 = 1.5 and 30; the top quark 
and the common squark masses are taken to be = 174 GeV and Ms = 1 TeV, 
respectively. The dashed line corresponds to the case where At = Ab = fi = (only 
the leading radiative correction is included), while the full lines correspond to the 
case where At = Af, = 1 TeV and fi = — 200, 0, +200 GeV (from top to bottom). 
The masses of the AiSSAi Higgs bosons h^, and H^, as a function of the 
pseudoscalar Higgs mass for the two previous values of tg/3 and for At = Af, = fj, = 
0,Ms = 1 TeV, are displayed in (b), (c) and (d), respectively. 

Fig. 15: The coupling parameters of the AiSSAi neutral Higgs bosons as functions of the 
pseudoscalar A° Higgs mass for two values of tg/? = 1.5 and 30 and for At = Aj, = 
/i = 0, Ms = 1 TeV. The couplings are normahzed to the SA4 couplings as defined 
in Table |1|. 

Fig. 16: The branching ratios of the AiSSAi Higgs bosons /i^ (a), (b), ^4° (c) and 

(d) as functions of their masses for two values of tg/3 = 1.5 and 30; the values 
At = Ab = fi = and Ms = 1 TeV have been chosen. [The arrows in (a) denote the 
branching ratios in the SAi limit of large A^ masses.] 

Fig. 17: The total decay widths of the AiSSAi Higgs bosons h^, H^, A^ and if^ as functions 
of their masses for two values of tg[3 = 1.5 (a) and tg/? = 30 (b); the values 
At = Ab = fi = and Ms = 1 TeV have been chosen. 

Fig. 18: Real and imaginary parts of the QCD radiative correction factor to the quark am- 
plitudes of the two-photon couplings for the AiSSAi neutral Higgs bosons: (a) h'^ 
and and (b) A'^; the renormalization scale for the quark mass is taken to be 
fiQ = m$/2. 

Fig. 19: The QCD corrected partial decay widths into two photons of the AiSSAi Higgs 
bosons h^, A^ for (a) tg/3 = 1.5 and (b) tgP = 30, and the the size of the 
QCD radiative corrections to the processes /H^ —>■ 77 and A'^ —> 77 (in %) as 
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functions of the Higgs boson masses for two values of tg/3 = 1.5 (c) and 30 (d). The 
renormahzation scale for the quark mass is taken to be — m$/2. 

Fig. 20: The deviation AE'-^ (a) and A£'^ (b) of the coefficients of the radiative QCD 
correction factors to the process ^ ^ gg from their values in the heavy quark limit, 
for two values of tg/? = 1.5 and 30; the renormahzation scale is taken tohe /i — m$. 

Fig. 21: The QCD corrected gluonic partial decay widths of the M.SSM. neutral Higgs 

bosons hP,H^ (a) and (b), for two values of tg/9 = 1.5 and 30; the size of 
the QCD radiative correction factor for h^/H^ gg (c) and gg (d). The 

renormahzation scale is taken to be = m$. 

Fig. 22: K factors of the QCD corrected cross sections (T{pp /H^ + -^) (a) and (j{pp 

A^ + X) (b) for two values of igf3 = 1.5 and 30; K^irt and K^b {A, B = q,g) arc the 
regularized virtual correction and real correction factors, respectively, and Ktot is 
the ratio of the QCD corrected total cross section to the lowest order cross section. 
The renormahzation and factorization scales are taken to he /i — M — m$ and the 
GRV parameterization for the parton densities have been used. 

Fig. 23: The dependence of the total K factors for the processes a{pp ^ + X) on the value 
of tgP for a characteristic set of Higgs boson masses. 

Fig. 24: The spread of the A4SSA4 Higgs production cross sections a{pp — > / + X) (a) 
and (t{pp A^ -\- X) (b) for two parameterizations of the parton densities. 

Fig. 25: The total production cross sections of the scalar CP-cvcn Higgs bosons h^,H^ (a) 
and the pseudoscalar Higgs boson A^ (b) at the LHC for two different cm. energy 
values: ^/s ^ U TeV and ^/i = 10 TeV. 

Fig. 26: The renormalization/factorization scale dependence of the M.SSM. Higgs boson 
production cross sections at lowest and next-to-leading order, for a characteristic 
set of Higgs boson masses and tg/3 values. 
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